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O . Abstract 

cn 

, The possibility for occurrence of a Bose condensate of dibaryons in nuclear mat- 

. ter is investigated within the framework of the Walecka model in the mean-field ap- 

0> ' proximation. Constraints for the w- and cj-meson coupling constants with dibaryons 

following from the requirement of stability of dibaryon matter against compression 
are derived and the effect of a- and 7r-meson exchange current contributions to the 

o ; 

^ , cr-dibaryon coupling constant is discussed. The mean-field solutions of the model are 

constructed. The effective nucleon mass m*^ vanishes when the density of dibaryons 
^ ■ approaches a critical value Pfyy^^ ~ 0.15 fm^^. The Green's functions of the equi- 

' librium binary mixture of nucleons and dibaryons are constructed by solving the 

Gorkov-Dyson system of equations in the no-loop approximation. We find that 
when the square of the sound velocity is positive, the dispersion laws for all ele- 
mentary excitations of the system are real functions. This indicates stability of the 
ground state of the heterophase nucleon-dibaryon mixture. In the model consid- 
ered, production of dibaryons becomes energetically favorable at higher densities as 
compared to estimates based on a model of non-interacting nucleons and dibaryons. 
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1 INTRODUCTION 

In 1977 R.Jaffe predicted the existence of the loosely bound dihyperon, H, with a 
mass just below the AA threshold. The H particle, if it exists decays only weakly. The 
exciting prospect to observe the long-lived H particle stimulated considerable theoretical 
and experimental activity. Calculations of the -ff-particle mass in the bag models 0, 0, 
the constituent quark model ^, the Skyrme model |0, on the lattice and in other 
models [§] showed that the existence of a dihyperon near the AA-threshold is plausible. 
The searches were proposed to examine the H particle production in proton-proton |^ 
p!0[] , proton-nucleus [|lT], |12|, nucleus-nucleus |jl2[ collisions, via a {K~,K^) reaction on a 



nuclear target |T^, and through a strong decay of a S"-atom system Weak decays 
of the H have also been studied [|r^ . The experiments |^, |T^, [T^ did not give a positive 
sign for existence of the H particle, however, a weak decay of the H produced in the 
p — C reaction has been reported [T^. The existence of the H particle remains an open 
question which must eventually be settled by experiment. The candidates for double- 
lambda hypernucleus whose existence constrain the binding energy of the H have been 
observed [ll9| . 

The non-strange dibaryons with exotic quantum numbers, which have a small width 
due to zero coupling to the A^A^-channel, are promising candidates for experimental 
searches The lowest-lying isospin T = dibaryons with quantum numbers 

JP = l-(iPi), S'i^Fs), l+{^Si,^Di), 2+(3D2),etc. couple to the A^A^-channel They have 
large widths and are difficult to observe experimentally. On the other hand, dibaryons 
made up of exotic quark clusters, for example, a and a cluster with relative or- 
bital angular momentum L = 1 may have unusual quantum numbers T = 0, = 0~. 
The data on pion double charge exchange (DCE) reactions on nuclei [^-[Q exhibit a 



peculiar energy dependence at an incident total pion energy of 190 MeV, which can be 
interpreted as evidence for the existence of a narrow d' dibaryon with quantum num- 



bers T = 0, = 0~ and the total resonance energy of 2063 MeV. Recent experiments at 
TRIUMPF (Vancouver) and CELSIUS (Uppsala) seem to support the existence of the d' 
dibaryon plj] . 

A method for searching narrow, exotic dibaryon resonances in the double proton- 



proton bremsstrahlung reaction is discussed in Ref. and some indications for a 

(ii(1920) dibaryon in this reaction have recently been found p8| . 

Dibaryons can be formed in nuclear matter. The properties of nuclear matter with 
admixture of multiquark clusters are discussed by Baldin et al. and Chizhov et al. [3U 



In these papers, the interaction between nucleons and multiquark clusters is included 
through a van der Waals volume correction. The occurrence of a heterophase state of 
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nucleons and 6-quark clusters is found to be energetically favorable in a wide region of 
temperatures and densities. A model for nuclear matter with an admixture of dibaryons 
with the short-range nuclear forces approximated by a 5-function-like pseudopotential 



is discussed by Mrowczynski |]3T[. Nuclear matter with a Bose condensate of dibaryons 
belongs to a class of heterophase substances whose properties are reviewed by Shumovskii 



and Yukalov 132 



The occurrence of a dibaryon Bose condensate in nuclear matter results in a softening 
of the equation of state (EOS) of nuclear matter. In the ideal gas approximation, the 
incompressibility of the heterophase nucleon-dibaryon matter vanishes. Occurrence of a 
dibaryon Bose condensate in interiors of neutron stars decreases the maximum masses 
of neutron stars P^l. The effect of the strongly interacting if-particle on the structure 



of massive neutron stars is investigated by Tamagaki |3J] and Olinto et al. The 
if-particle interaction with nucleons and the ifif-intercaction are studied in the non- 
relativistic quark cluster model 0, . 



In a recent paper an exactly solvable model for a one-dimensional Fermi-system of 
fermions interacting through a potential leading to a resonance in the two-fermion channel 
is constructed. This model takes the Pauli principle for fermions and the composite nature 
of the two-fermion resonances into account. The behavior of the system with increasing 
density can be interpreted in terms of a Bose condensation of two-fermion resonances. 

A Bose condensate of dibaryons can presumably be formed at higher densities when 
relativistic effects become important. In order to describe such a system, one should 
go beyond the non-relativistic many-body framework. The relativistic field-theoretical 
Walecka model ^ is known to be very successful in describing properties of infinite 
nuclear matter and of ordinary nuclei throughout the periodic table. It constitutes the 
basis for the quantum hadrodynamics (QHD) approach for studying nuclear phenomena. 

In this paper we study the effect of narrow dibaryon resonances on nuclear matter in 
the framework of the Walecka model. The Lagrangian of the model contains nucleons 
interacting through uj- and a -meson exchanges. We add to the Lagrangian dibaryons 
interacting with nucleons and with each other via the exchange of uo- and a-mesons. 

The outline of the paper is as follows. In the next Sect., we start with a discussion 
of the properties of heterophase nucleon-dibaryon matter in the ideal gas approximation. 
We give simple qualitative estimates at what densities (chemical potentials) and for which 
dibaryon masses a dibaryon Bose condensate may occur. The softening of the EOS for 
neutron matter and its consequences for the gravitational stability of neutron stars is 
discussed. In Sect. 3 an extension of the Walecka model including dibaryon fields in 
the Lagrangian is presented. We discuss in detail the constraints for the meson-dibaryon 
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coupling constants which follow from the requirement of stability of dibaryon matter 
against compression. We also calculate the numerical values for these coupling constants 
in the additive model and estimate the size of o- and tt- meson exchange current corrections 
to the cr-dibaryon coupling constant. 

In Sect. 4, the mean-field (MF) solutions to the extended Walecka model are con- 
structed. We show that the self-consistency equation for the effective nucleon mass in 
presence of the dibaryon component can formally be reduced to an analogous equation of 
the standard Walecka model. The MF solutions are examined from the point of view of the 
Hugenholtz-Van Hove theorem |Q which serves as a check for the internal consistency. It 
is shown that the thermodynamic pressure coincides with the hydrostatic pressure. The 
numerical results show a softening of the EOS for nuclear matter with a dibaryon com- 
ponent. In Sect. 5, we construct Green's functions of the heterophase nucleon-dibaryon 
matter by solving a system of the Gorkov-Dyson equations in the no-loop approximation. 
We find that when square of the sound velocity is positive, the dispersion laws for ele- 
mentary excitations are real functions, indicating the stability of the heterophase ground 
state of nucleon-dibaryon matter. 

2 IDEAL GAS APPROXIMATION 

In the ideal gas approximation, the physical picture of dibaryon condensation is very 
simple. Suppose we fill a box with neutrons (see Fig.l). Due to the Pauli principle, the 
neutrons occupy successively higher energy levels. This process is continued until the 
chemical potential of two neutrons on top of the Fermi sphere is larger than the dibaryon 
mass. When the Fermi energy of the neutrons becomes larger, it becomes energetically 
favorable for the neutrons to form dibaryons. The critical density for dibaryon formation 
is determined by the mass of the lightest dibaryon. Above the critical density, the chem- 
ical potential of the nucleons /i„ is frozen at the value /x™^"^ = m£)/2 where is the 
dibaryon mass. This equation is a consequence of the chemical equilibrium with respect to 
transitions nn D. Dibaryons are Bose particles, so they are accumulated in the ground 
state and form a Bose condensate. Dibaryons have zero velocities, therefore they do not 
collide with the boundary and do not contribute to the pressure. Because the Fermi en- 
ergy of the neutrons is frozen, the pressure does not increase with the density (see Fig. 
2). Consequently, nuclear matter loses its elasticity and its incompressibility vanishes. 
Nuclear matter with such properties cannot protect neutron stars against gravitational 
compression and subsequent collapse. 

Obviously, a dibaryon Bose condensate does not exist in ordinary nuclei. The equality 
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fJ-n = f^p follows from the charge symmetry of nuclei. Assuming that the shell model 
potential for dibaryons is twice as deep as the one for nucleons, we conclude that the 
masses of dibaryons coupled to the iViV-channel should be greater than 

mo > 2fiN = 2(m7v + Sp) = 1-96 GeV (2.1) 

where ep = 40 MeV is the Fermi energy of nucleons in nuclei. For example, the d' dibaryon 
[^-[^ is coupled to the NNn channel only. In the medium, the reaction pd' ^ nnp 
is possible. In nuclei the equilibrium condition for the chemical potential has the form 
fiD + f^p = 2/i„ + /ip. Since /ip = /U„ = ^n, we arrive at the same constraint (^TTP . 

In chemical equilibrium with respect to the /3-decays p + e ^ n + u^, the chemical 
potentials for neutrons and protons in neutron matter satisfy the equality: yU„ = fip+fie, so 
that dibaryons coupled to the nn-, np-, and pp-channels actually occur at slightly different 
densities; and dibaryon condensation starts if dibaryon masses are smaller, respectively, 
than 2yU„, 2yU„ — yUg, and 2yU„ — 2/ie. In case of the d', the critical value for the neutron 
chemical potential is determined by the condition m^/ = 2/i„ — fie- The electron chemical 
potential fie is positive. The H particle is coupled through a double weak decay to the 
A^A^ channel so that the H particle occurs at rriH = 2/i„. 

Nonrelativistically, the density p and pressure p distributions inside of neutron stars 
are described by Euler's equation 

dv 

p— + p(v ■ V)v = -Vp - pV$ 
ot 

where $ is the gravitational potential and v is the velocity of the neutrons. If a dibaryon 
Bose condensate is formed, the pressure in the internal region of the neutron star should 
be constant (the dotted horizontal line on Fig. 2), and therefore Vp = 0. Gauss's law 
implies / dS ■ V$ = 47rGM(r), where M(r) is the mass inside of a sphere of radius r. It 
follows that V$ 7^ 0. From Euler's equation, we get v 7^ 0, i.e. there is no static solution. 
We thus conclude that in the ideal gas approximation there are no stable solutions if 
a dibaryon Bose condensate is formed and neutron stars are gravitationally unstable. 
These physical arguments are also valid in general relativity where the neutron stars are 
described by the Oppenheimer-Volkoff equation and are qualitatively correct beyond 
the ideal gas approximation Therefore there is an interesting connection between 

the masses of the lightest dibaryons and the upper limit for masses of neutron stars. 

The masses of several neutron stars are reliably determined to be above IAMq. In 
the tensor interaction model (a stiff model) and in the Reid model (a soft model; for a 
review of these nuclear matter models see e.g. |^5[)' neutron chemical potential in the 
center of a mass I.AMq neutron star can be evaluated to be 1090 MeV and 1015 MeV, 
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respectively. The requirement that there be no Bose condensate of dibaryons coupled to 
the iViV-channel (like the H -particle) inside neutron stars gives 



mo > 2.18 GeV and tud > 2.03 GeV. (2.2) 

These numbers are in the range of present experimental searches for dibaryons. They are 
valid, however, only when the interaction of dibaryons is neglected. 

The interaction of dibaryons with neutrons and with each other increases the pressure. 
The equation of state of neutron matter becomes stiffer, yielding stability of neutron stars 
in some interval of densities [^]. If dibaryons are formed in a first order phase transition 
(such a scenario is discussed by Tamagaki |Q for the if-particles) , neutron stars become 



unstable at a critical density provided that the jump Ap of the density in the phase 
transition point is sufficiently large ^]. In Fig. 2 Ap = oo. In Newtonian gravity, 
the criterion is given by Ap > (3/2)p. 



3 DIBARYON EXTENSION OF THE WALECKA 
MODEL 

Many successful phenomenological applications of QHD demonstrate that the interactions 
of hadrons at large and intermediate distances can be adequately described in terms of 
hadronic degrees of freedom. Many observables are not sensitive to the contributions from 
very short distances. In QHD, the effects of retardation and causality are rigorously taken 
into account. 

Parameters of the Walecka model are fixed by fitting the properties of nuclear matter 
at the saturation density. Once the model parameters are fixed, other consequences can 
be extracted without any further assumptions. The inclusion of dibaryons to the model 
entails several uncertainties due to the lack of reliable information on dibaryon masses and 
coupling constants. However, many conclusions can be drawn on quite general grounds 
without knowing precise values for the newly added parameters. 



3.1 Lagrangian of the model 

We consider an extension of the Walecka model by including dibaryon fields to the La- 
grangian density 

L = ^{id^'y^ -ruN- 9^(7 - ^..c^m7m)^ + U^t^cr)^ - \mla'^ - \Fl^ + |m>2+ 
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Here, \l/ is the nucleon field, uj^ and a are the u- and a-meson fields, F^^ = d^uj^ — d^Uy 
is the field strength tensor of the vector field; is the dibaryon field, for which we assume 
that it is a isoscalar-scalar (or isoscalar-pseudoscalar) field. This assumption includes the 
interesting cases of the if -particle and the d'-dibaryon. The values and mo- are the uj- 
and cr-meson masses and Qu)-, Qa-, h^j, h^r are the coupling constants of the u- and a-mesons 
with nucleons (g) and dibaryons (h). 

The term Cc describes the conversion of dibaryons into nucleons. The ii-particle is 
coupled to the A^A^-channel through a double weak decay, so that Cc = 0{G'j;^). For 
the non-strange di and the d' dibaryon, we neglect possible virtual transitions e.g. to 
the NNa channel. The on-shell couplings for these dibaryons are small. The exotic di 
dibaryon decays to the A^A^7-channel only, and therefore Cc = 0{a). The d' dibaryon 
decays to the NNn channel. Due to the Adler's self-consistency condition Cc oc d^ir. 
In the MF approximation S^vr = 0, and the term Cc does not modify the MF equations. 
In what follows we set Cc = 0. The effect of a small term Cc reduces to providing a 
chemical equilibrium with respect to transitions between dibaryons and nucleons. 

The field equations corresponding to the Lagrangian have the form 

{id^.'Jf, -rriN- ga<y - S-c^t^^T/.)^ = 0, (3.2) 

(-□ - ml)a = g^^^ + 2K{mD + Ko)^*^, (3.3) 

((-□ - ml)g^,^ - d^,dy)ujy = + h^(p*i ^^^(p- 2hlu^(p*ip, (3.4) 

((•9^ + ih^oOf^y + {niD + K(yf)<^ = 0. (3.5) 
The field operators can be expanded into c-number- and operator parts: 

a = ac + a, 

(3-6) 

Lp = LPc + V^ 

The c-number parts of the fields A = a, u^, (p, and ip* are defined as expectation values 
Ac =< A > over the ground state of the system. The average values of the operator parts 
are zero by definition: < A >= 0. 

The cr-meson mean field determines the effective nucleon and dibaryon masses in the 
medium 

m^f = niN + ga(^c, (3.7) 
m*^ = niE, + h^ac- (3.8) 
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The nucleon and dibaryon currents have the form 

= ^7m^, (3.9) 

jj^ = ip*i ^^l - 2h^LJ^ip*(p. (3.10) 
The baryon number current is defined by 

^ =^ + 2^. (3.11) 

The cj-field is coupled to the current 

= 9.3^, + (3.12) 
The nucleon vector and scalar densities are defined by expectation values 

Pnv =< ^7o* >, (3.13) 

Pns =< > . (3.14) 
The scalar density of the dibaryon condensate is defined by 

Pd5=I<¥'(0)>|'. (3.15) 

The time evolution of the condensate 99- field is determined by the chemical potential po 
of dibaryons 



V,{t) = e-^^-'^p%s- (3.16) 

It is useful to separate the contribution of the a;-meson mean field to the chemical 
potential energy of dibaryons 

Pd = P*D + hujUJc- (3.17) 
The dibaryon number density is according to Eg. (|3.10|) given by 



p'hv = 2p},p^j,s- (3.18) 

The possibility of existence of a Bose condensate of dibaryons depends on values of the 
coupling constants of dibaryons with the u- and cr-mesons. 

3.2 Stability of dibaryon matter against compression 

The uj- and a- meson coupling constants h^j and enter the Yukawa potential for two 
dibaryons 

V{r) = . (3.19) 

47r r An r 
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The interaction energy for dibaryons in the condensate for a constant density distribution 
P£)(x) = pd is equal to 

W=^l rfxirfx2PD(xi)pi5(x2)K(|xi - X2I) = 27riVz)PD(^^ - (3.20) 

where Nd is the total number of dibaryons. The integral (3.20) is linear in density. A 
negative W would imply instability of the system against compression. The value W is 
positive and the system is stable for 

h"". hi , , 

> 7-±7- (3-21) 



In a nonrelativistic theory for interacting bosons and in the model considered (see 



Sect. 4), the requirement of stability is equivalent to the requirement of a negative value 
for the boson forward scattering amplitude and/or a positive value of square of the sound 
velocity (a^ > 0). In the Born approximation, the forward scattering amplitude is, as in 
Eq.(3.20), proportional to the volume integral of the potential. 

The iZ-particle interactions were studied in the non-relativistic quark cluster model 
0, |26[-|^ which is successful in describing the A^A^-phase shifts. The calculation of the 



interaction integral ( |3.20D with the adiabatic if if -potential gives a negative number, 
so that the ii-dibaryon matter is probably unstable against compression. The coupling 
constants of the mesons with the ii-particle can be fixed by fitting the depth and the posi- 
tion of the minimum of the iiii-potential to give /i^ = 603.7 and /i^ = 279.2. These values, 
while yielding a negative W ^ probably overestimate the repulsion between ii-particles at 
small distances, and therefore underestimate \W\. The meson-dibaryon coupling constants 
for the case of the d\ and d dibaryons are presently unknown. 

3.3 Coupling constants in the additive model 

In the additive picture, mesons interact with the constituents of the dibaryon (Fig. 3). 
For non-strange dibaryons coupled to the A^A^-channel, the a - and uj- meson couplings 
are in the nonrelativistic approximation simply twice the corresponding meson-nucleon 
coupling constants: = 2g^ and = 2ga. The scalar charge is, however, suppressed by 
the Lorentz factor. This effect decreases the value of h„. Note that the magnitudes of the 
meson coupling constants for the if -particle, extracted from the adiabatic if if -potential 



are consistent with the additive estimates: hi^/{2g^) = 0.89 and h^/{2g^) = 0.80. 



For the standard set of parameters of the Walecka model [|T, rricr = 520 MeV, 
gl = 190.4, and gl = 109.6, the inequality (^72T\ ) 

98.85(-^)2GeV-2 > i29.0(-^)2GeV~2 (3 22) 

Sfi-o; 2g^ 
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is not satisfied. The precision of the additive estimates is, however, not better than 30% 
of the central values. Exchange current contributions which violate the additivity are 
discussed below. 



The violation of the inequality ( 3.21 ) for the additive estimates of the meson-dibaryon 
couplings is not accidental. For homophase nuclear matter, g^^ should be greater than 
Qa in order to get sufficient repulsion between nucleons at small distances. To reproduce 
the properties of nuclear matter at the saturation density, the following inequality should 
hold 

<7^- (3-23) 



At small densities, the inequality p.23| ) provides a negative value for the interaction in- 
tegral (3.20) between the nucleons, resulting to the local instability of nuclear matter 
against compression. When the density increases, the relativistic effects become impor- 
tant. The scalar density of nucleons increases slower than the vector density, since it is 
suppressed by the Lorentz factor < I/7 >. This finally leads to an equilibrium at the 
saturation density of nuclear matter. 

For exotic dibaryons like the rf', one should take into account the presence of a pion 
in the resonance wave function. The cu-meson decays into three pions. It is not coupled 
to the pion in the d! wave function (see Fig. 3), so = 2g^. 

The 0" — TT cubic couplings are described by the effective Lagrangian density 

AC = --a^--an^ (3.24) 
6 6 

where tt is the pion field. The cubic terms generate three-body forces between nucleons. 
Phenomenological fits to the bulk nuclear properties give 



0.9^5.3. (3.25) 

rriN 

The lower and upper values correspond, respectively, to a small negative and a large 
positive term Xa^ in the effective Lagrangian. 

The non-linearities of the linear sigma-model are qualitatively different [25] 

K k' mt — mt 
= = -3^^ " „ " = -8.9. (3.26) 

In the additive model for d', we get 

= 1 + = 1 + 0.027—. (3.27) 

2g„ 2Ag„m^ mN 

The last term describes the contribution from the first of the diagrams in Fig. 3(b). For 
\k' /mN\ < 10 the correction to the additive value is smaller than 30%. The sign, however, 
is not defined. 
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The Brown-Rho scaling for non-strange hadron masses is reproduced at the tree 



level for = {mD/mN)ga- In such a case, 



rriisi mj:, 

Since niE, > 2m at, Brown-Rho scaling gives = {Tn£,/m]^)ga^ > 2g„. 



(3.28) 



3.4 Exchange current contributions to the coupUng constants 

The exchange current contributions to the a-dibaryon coupling constants shown in Fig.4 
can be extracted from the Lagrangian (3.24). The a-meson field generated by a pointlike 
source of charge h„ is determined from the equation 

(A-m^)a(x) = M(x). (3.29) 

The sign of the right hand side of the equation is chosen such as to yield a negative 
(attractive) a-meson field for a positive h„. It is useful to compare Eq.(3.29) with the 
static limit of the equation of motion for the cx-meson field determined by the Lagrangian 
C + AC with no dibaryon component: 

(A - m^)a(x) =g^Y: ^'fc(x)M/,(x) + ^^^(x) + '^n^^). (3.30) 
fc=i " 

The constant measures the scalar charge generating the a-meson field around dibaryons, 

so that we can write 

2 / 

K = [ d^[g^ J2 ^fc(x)^.(x) + ^a2(x) + ^7r2(x)]. (3.31) 
^ k=i 2 6 

Nonrelativistically, the cr-meson field created by nucleons located at points x/^ is 

2 g-mCT|x-Xfc| 

a{^) = -J2g^— -. (3.32) 

47r|x-Xfc| 

Substituting this expression into Eq.(3.31) and omitting diagonal terms, we get for the 
exchange-current contribution an expression 

AK{Dy~^^^ = -^^gl < e-^-l'i-'^l > (3.33) 

where xi and X2 are coordinates of the two nucleons. The exponential term should be 
averaged over the dibaryon wave function. 

To give an order-of- magnitude estimate, we use for the relative wave function of two 
nucleons Hulthen wave functions 
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with /i = 1 f m ^ and v = 4/i. The correction (3.33) can be estimated to be 

— = 0.040 . 

2g^ iriN 

The d'-dibaryon contains a constituent 7r-meson. The a-meson field created by the 
pion is given by 

a(x) = — (3.35) 

^ ^ 12m^47r|x-X3| ^ ^ 

where X3 is the pion coordinate. The total change in the a -meson coupling with the d' 

equals 



lovrmo- 



6m^ 



(3.36) 



Assuming that the two-body probability densities are given by the square of the wave 
function (3.34), we get 

= 0.040 1 + 0.108 . 

The contribution of the pion exchange currents can be evaluated in a similar way. The 
pion field created by nucleons has the form 



2 



^(x) = E ^T^Ma, ■ -. (3.37) 

t^i 2mN 47r|x-Xfc| 

For dibaryons consisting of two nucleons and also for the d', the exchange pion current 
contributions are given by 

1447r 2mN \y|xi— X2I / / 

(3.38) 

The evaluation of this expression using the wave function (3.34) gives 

— ^LJ = 0.0015 (n ■ r2)(ai ■ as). 

2ga rriN 

The pion exchange current correction to the a-dibaryon coupling is therefore significantly 
smaller than the corresponding sigma exchange current corrections. 

In view of the estimates (3.25) and (3.26) the exchange current corrections cannot 
exceed 30% — 40% of the additive values. At present, the coupling constants of the mesons 
with dibaryons are not known with sufficient precision to draw any definite conclusion 
concerning the stability of dibaryon matter. 

The if-particle, however, is studied better than other dibaryons. In what follows, we 
use a realistic if if- interaction, based on the quark cluster model p6|-|l39|. 



g-m^|x-Xfe 



12 



4 THE MEAN-FIELD SOLUTIONS 



The MF solutions of the Walecka model are asymptotically exact in the high density limit 
[^ . They serve as a starting point for the calculation of corrections for finite densities 
within the QHD approach. 

Neglecting the operator parts of the meson fields in Eqs.(3.2)-(3.5), we get the following 
expressions for the meson mean fields 

mi mi 

_ gapNS + K2mDp'i)s _ g^pNS + K2m})p'j^s 
m^ m^ 
The effective masses of mesons in a heterophase nucleon-dibaryon matter are given by 

ml = ml + 2hlp'hs, (4.3) 

ml = ml + 2hlp'ts. (4.4) 

The mechanism responsible for the change of meson masses is essentially the Higgs mech- 
anism. 

Substituting expression (3.16) into Eq.(3.5), we get 

P*D = (4.5) 

The nucleon and dibaryon chemical potentials have the form 

UN = E*p + 9u^c, (4.6) 

PD = m}) + KuJc, (4.7) 



where Ep = dm*^ + kp is the Fermi energy of nucleons with the effective mass m^. 



4.1 The self-consistency equation 

The nucleon vector and scalar densities are given by 

r d\Q 

PNV=1 77r^0ikp-\k\), (A.i 



. (27r)3 



The statistical factor 7 = 4 (2) for nuclear (neutron) matter. 
The self-consistency equation has the form 



m*N = mN - ^{gaPNS + Kp'hv)- (4-10) 
mt 
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It follows from Eqs.(3.7), (3.18), (4.2), and (4.5). The total baryon number density 
equals pxv = Pnv + '^Pdv- Eq-(4.10) can be transformed to a form equivalent to the 
self-consistency equation of the standard Walecka model without dibaryons: 

9^ 

m*]^ = mN ^Pns, (4.11) 

where 

c,max c 

m;v = m^ ^^^,,~f^^ , (4.12) 
Pdv 

pi)y = — = 0.1507(— )fm . (4.13) 

If the densities ptv and p^y are fixed, equation (4.11) allows to find the effective nucleon 
mass m^. Solutions to Eq.(4.11) exist for an arbitrary total density pxv when the value 
ifiN \s positive. This is the case for p^y < p'}^^^. Note that the dibaryon mass does not 
enter Eq.(4.11) directly. 

In Figs. 5 (a,b) we show for = 2g^^ the critical density for occurrence of a Bose 
condensate of dibaryons in (a) symmetric nuclear matter and (b) neutron matter as a 
function of the cr-dibaryon coupling constant. The critical density is determined from 
equation 2p7v — Pd = 0. In the region > the requirement (3.21) is fulfilled. The 
dibaryon components of the heterophase nuclear- and neutron-dibaryon matter are stable 
against compression. 

For 2mN <i Tna < 1.89 GeV, we start at zero density from heterophase nuclear- 
dibaryon matter (a). With increasing density, the matter transforms into homophase 
nuclear matter and then again back to heterophase nuclear-dibaryon matter. For rriD > 
1.89 GeV, we start at zero density from homophase nuclear matter which converts with 
increasing density (at prv > Po for ha/{2ga) < 0.8754) into heterophase nuclear-dibaryon 
matter. 

When the cr-meson coupling with dibaryons decreases, the effective dibaryon mass 
m|) increases (see Eq.(3.8)) and dibaryon formation is therefore suppressed. We thus 
conclude that the curves in Figs. 5 (a,b) should have a negative slope for a transition from 
homophase to heterophase matter and a positive slope for a transition from heterophase 
to homophase matter. The region of instability for neutron matter is greater because at 
the same density we have higher values of the neutron chemical potentials and therefore 
more favorable conditions for the production of dibaryons. 

We show the results for rfi-dibaryon |^ with possible quantum numbers = 1"*". For 
such dibaryons, the MF equations derived above for the scalar (pseudoscalar) dibaryons 
jp = 0^ remain the same. In case of the = 1^ dibaryons, the dibaryon part of the 
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Lagrangian ( |3.1| ) can be written in the form 

Ho = -D^iflD^if, + {mo + KaYif^if^ + XD^^l + \*D^ip, 
where = + ihujUjfj, and A and A* are Lagrange multiphers. The field equations are 

{Dl + {mn + Kaf)^^ - D^X = 0, 

The last equation removes the unphysical timelike component of the dibaryon vector field 
in the co-moving frame of the particle. It is equivalent to the requirement k^v^^ = where 
Un is a four-vector velocity of the particle. 
The condensate field evolves like 



where e is a unit vector. The field equations give A = A* = and result in Eqs.(4.1) and 
(4.2). The energy and pressure have the form (4.22) and (4.23) of the = 0^ dibaryons. 

In Fig. 6 we classify possible behaviors of the difference 2fiN — fi^ between the chemical 
potentials of nucleons and dibaryons with growth of the dibaryon fraction. 

It is clear that when the difference is positive and the dibaryon density is zero, p^y = 
0, production of dibaryons is energetically favorable. In such a case, the fraction of 
dibaryons increases. The state pjjy = is therefore unstable. If the difference 2/iAr — 

is negative and the substance consists of dibaryons only, production of nucleons is 
energetically favorable. This state is unstable. If the difference 2fi^ — p^ is zero, but 
increases with the dibaryon fraction, small fluctuations bring the substance away from 
the equilibrium. Such a state is unstable also. 

The system is stable in the following three cases. 

(i) Homophase nuclear and neutron matter: 

(ii) Homophase dibaryon matter: 

(4-15) 

2pz)y = Ptv- 

(iii) Heterophase nucleon-dibaryon matter: 

- /^D = 0, (4.16) 

d{2pN - Pd) I 



dpDV 



\PTV 



< 0. (4.17) 
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In the first case there are no dibaryons, in the second case there are no nucleons, and in 
the third case we have a heterophase mixture of nucleons and dibaryons. Small fluctua- 
tions around the state 2fi]\f — fi^ = bring the substance back to the equilibrium point. 
Eqs.(4.16) and (4.17) therefore describe a stable equilibrium. 

In Fig. 7 (a) we show the nucleon effective mass m*j^ versus the dibaryon fraction 
'^Pdv/Ptv in heterophase nuclear matter and in heterophase neutron-dibaryon matter for 
the coupling constants /i^ = 2(7^ and h(j/{2g„) = 0.8. At the same total baryon number 
density, neutron matter is more relativistic. The scalar charge density of neutron matter 
is thus lower, the ci-meson mean field is smaller, and therefore the effective nucleon mass 
is greater. On the plot, the dotted lines corresponding to neutron matter lie above the 
solid lines corresponding to symmetric nuclear matter. 

When the dibaryon vector density p^y approaches its maximum value p^y^"", the 
effective nucleon mass vanishes. This effect can be interpreted as follows. Two nucleons 
on the top of the Fermi sphere have energy 2Ep. In chemical equilibrium with respect to 
transitions NN <-> D, the relation 2Ep — holds. In transitions NN <-> D the baryon 
vector charge does not change. However, the scalar charge changes. For two nucleons 
the scalar charge equals 2g^^, whereas a dibaryon in the condensate has scalar charge 
ha-- When the system is nonrelativistic, formation of new dibaryons is accompanied by a 
decrease of the scalar charge density, since 20^^ ~ 2ga > I-GQct — h^. This phenomenon 
is reflected in the slight increase of the effective nucleon mass with the dibaryon fraction 
at low total baryon number densities. When the density is high, the system becomes 
relativistic, and so the Lorentz factor comes into play. The scalar charge of the two 

^F 

nucleons is small, whereas the dibaryon scalar charge is large. As a result, the scalar charge 
density increases, the scalar mean field increases, and the effective nucleon mass decreases 
with the dibaryon fraction. In the standard Walecka model, the effective nucleon mass 
vanishes at infinite density. In the model considered, the effective nucleon mass vanishes 
when the density of dibaryons approaches the value Pjy^^^ . Note that the behavior of the 
effective nucleon mass m*^ with growing dibaryon fraction, shown in Fig. 7 (a), does not 
depend on the vacuum value of the dibaryon mass, because the dibaryon mass does not 
enter the self-consistency equation (4.11). 

In Fig. 7 (b) we show the difference for the chemical potentials versus the dibaryon 
fraction for itld = 1.96 GeV. Since the vacuum dibaryon mass does not enter the self- 
consistency equation (4.11) and enters linearly in the difference 2fiN — P-d (see Eqs.(3.8) 
and (4. 6), (4. 7)), the results for other dibaryon masses can be obtained simply by vertical 
parallel displacements of the curves. The results for the mo — 2.08 GeV {d' dibaryon) 
can be obtained e.g. by a 100 MeV negative shift, etc. 
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The MF solutions exist at all densities prv for sufficiently small densities of dibaryons, 
P'dv < P^Dv^- K means that we can always investigate the stability of homophase nuclear 
and neutron matter with respect to formation of a dibaryon Bose condensate. When the 

total density Ptv is very high, the dibaryon production is energetically favorable. The 
MF solutions disappear, however, before the system reaches equilibrium. 

4.2 Thermodynamic consistency checks for the mean-field solu- 
tions 

The canonical energy-momentum tensor corresponding to the Lagrangian density (2.1) 
can be written in the form 

rp rpN 1 rpa 1 rpuj i rpD 

-I- fJ-i^ — -'- fii^ ^ -'^ HU ' ^ ' -'■ 1X1/ 



where 



rf, = m^,d,^, (4.18) 



1 



T^u = -^9i.M<^dra - m'y) + d^ad^a, (4.19) 

T^^u = ^gnAidr^^x - dxUJr)drUJx - mlurUr) - d^u^d^Ur (4.20) 
T^^, - 2d^^*d,^ - Kuj^^H du V. (4.21) 
The energy density e —< Tqq > given by average value of the Tqq component has the form 

£ = 7/^ j^{E*{k) + g^uc) + {m*j, + /i^u;c)p^y + i^mlal - -mlul (4.22) 

Here, E*{\s) + Quji^c is the total energy of a nucleon with momentum k, and + h^^uOc 
is the total energy of a dibaryon in the ground state (with zero total momentum) in the 
external a;-meson mean field. The last two terms are the contributions of the classical u- 
and (T-meson fields to the energy density. 

The hydrostatic pressure p = — | < Tjj > has the form 

Because dibaryons in the condensate are at rest they do not contribute to the pressure. 
In agreement with the general requirements 
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The pressure can be calculated in two different ways: from Eq.(4.23) and from the ther- 
modynamic relation 

p = -e + ixnPnv + fJ'DpDV (4-26) 

It is not difficult to check that the thermodynamic pressure (4.26) coincides with the 
hydrostatic pressure (4.23). 

The Hugenholtz-Van Hove theorem |Q (HV) requires that the energy of fermions at 
the Fermi surface and the average energy of a physical system at zero pressure (at the 
saturation density) are equal. This theorem is useful for checking the internal consistency 
of approximations. The mean-field theory and the relativistic Hartree approximation of 
the standard Walecka model obey this theorem. 

The vacuum dibaryon mass enters our model as a free parameter. The formal as- 
sumption that the dibaryon condensate occurs at the saturation density does not lead to 
any contradictions. In chemical equilibrium with respect to the reaction A^A^ ^ D, the 
relations (4.16) are fulfilled. In agreement with the HV theorem, for p = Eqs.(4.6) and 
(4.26) give 

fiN = E*p + g^Uc = (4.27) 
Ptv 

In Figs. 8 and 9 we show the energy per nucleon and the pressure versus the total 
baryon number density for some possible dibaryons. It is useful to compare Figs. 2 and 9. 
In Fig. 2, we schematically show the behavior of the pressure as a function of the density in 
the ideal gas approximation. In Fig.9, the pressure is calculated for the dibaryon extension 
of the Walecka model. The effect of vanishing incompressibility in the ideal gas model 
(the horizontal dotted line in Fig. 2) is displayed in Fig.9 as a softening of the EOS of the 
heterophase nuclear (dashed line in Fig.9) and neutron-dibaryon matter (dotted line in 
Fig.9). The occurrence of if-dibaryons provides a possible mechanism for the formation 
of strange matter. 

Note that the pressure of the heterophase system of nucleons and dibaryons obeys the 
basic inequality of statistical mechanics [ p^ 

dp 



dpTV 



> 0. (4.28) 



4.3 The concept of equilibrium of heterophase substances for 

jCc = ^ and jCc 7^ 0. 

We discuss here the effect of a small term Cc describing transitions between dibaryons 
and nucleons. 

In homophase substances, the energy density is a function of the total baryon number 
density prv only. 
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e = e{pTv)- 

In a heterophase substance, the energy is a function of additional parameters 



The parameters characterize components of the heterophase substance. The equihbrium 



state is determined by the conditions 

de{pTv,Ci) 
Eigenvalues of the matrix 



(A; = 1,2,...). (4.29) 

(4.30) 



should be positive to guarantee a minimum of the energy. 

In our case, the state is specified by the nucleon chemical potential pn, the dibaryon 
chemical potential pn, and by the expectation value of the dibaryon field p'j^g. The 
nucleon chemical potential p^ is a free parameter. The dibaryon chemical potential po 
determines the evolution of the condensate part of the dibaryon field (Eq.(3.16)). 

The value /id is fixed by the Hugenholtz-Pines condition [^. In the MF approxima- 
tion, p*j^ = m*jy. There remain two free parameters: pn and pj^g. The total density is a 
sum of the two terms prv = Pnv + '^Pdv- This expression is also valid for £c 7^ (if there 
are no derivatives in Cc). 

We can chose as independent parameters the values p^y and p^y. It is not necessary 
to require that pnv and p'j^y be timelike components of two conserved currents. It is 
sufficient that they characterize uniquely the phases of the binary mixture. Equations 
(|]29D give then 

^ de{pNv,phv) ^ de{pNv,p'Dv) 
dpNV dpDv 
One can use expression ([4.22|) for the energy density to verify that Eq.( [4.16|) follows 

from Eq.( 4.31 ). This is a consequence of Eqs.( [4.24 ) and ( 4.25 ). The matrix ( [4.30| ) is 



positively definite when Eq. (|4.17| ) is fulfilled. 

For every conserved current one can introduce an independent chemical potential. 
For £c = 0, the parameters pn and pd have the meaning of the chemical potentials 
corresponding to two conserved nucleon and dibaryon currents. For £c 7^ 0, there is only 
one conserved baryon current and only one baryon chemical potential. 

In the both cases (£c = and £c 7^ 0), the equilibrium is determined by Eq.( ^.31| ). 



For Cc = 0, this equation reduces to Eq.( |4.16 ). It can be naturally interpreted in terms 



of the chemical equilibrium between nucleon and dibaryon phases. 
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For Cc 7^ 0, the values p^qy and p^y are no longer timelike components of conserved 
currents. They play the role of formal parameters characterizing the nucleon and 
dibaryon phases. In such a case, the condition ( 4.31| ) cannot be interpreted in terms of a 
chemical equilibrium, since we cannot determine individual contributions of nucleons and 
dibaryons to the total baryon number of the system. 

For Lc 7^ 0, one should write additional terms 0(£c) on the right hand side of 
Eqs.(4.16) and (4.17). For narrow dibaryons which we discuss here, the corrections to 
these equations are small. For example, in case of the if-particle the corrections are of 
order 10~^°. 



5 GREEN'S FUNCTIONS AND ELEMENTARY EX- 
CITATIONS 

In the MF approximation, all dibaryons are in the Bose condensate. It is known that 
interaction between Bose particles brings some fraction of bosons out of the condensate. 
In the nonrelativistic approximation, the density of bosons which are not in the condensate 



increases with the total density faster than the density of bosons in the condensate |^ . 

The bosons that are not in the condensate have finite momenta and they contribute 
to the pressure. One expects that the dibaryons which are not in the condensate shift the 
critical density for disappearance of the effective nucleon mass to higher values, because 
the dibaryon scalar charge is suppressed by the Lorentz factor. The contribution of 
dibaryons out of the condensate to the energy density and pressure can be calculated 
using the diagram technique developed for Bose systems by Belyaev |5^. To go beyond 



the MF approximation, it is necessary to determine the Green's functions of the system. 
5.1 Solutions of the Gorkov-Dyson equations 

To eliminate the time dependence from the condensate parts of the (/^-fields, we pass to 
the //-representation for dibaryons. This can be done by the substitution (y9 — > (pe^^^^ and 
if* v9*e-^^^*. 

The Green's functions are defined by 

iG^p{x' - x) =< T^„(a;')^/3(x) >, (5.1) 

iD^^'ix' -x) = zD^^(x - x') =< tA{x')B{x) >, (5.2) 
with A, B = a, ui^, Lp* . In momentum space 

D^^(fc) = D^'^i-k). (5.3) 
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The (7- and a;-vertices with dibaryons are depicted in Fig. 10. The crosses on the 
dibaryon double hnes denote the appearance or disappearance of dibaryons from the Bose 
condensate. In Fig. 11 wc graphically show a representation of Eqs.(4.1) and (4.2). The 
set of these diagrams can be summed up either to modify the meson propagators (the 
bold dashed lines) or the meson vertices with dibaryons. The dressed meson-dibaryon 
vertices are determined by the diagrams shown in Fig. 12. 

It is useful to distinguish three kinds of propagators. The bare ones denoted by thin 
lines, the mean field ones denoted by thick lines, and the complete ones denoted by blobs 
with outgoing thick lines. The MF propagators are determined by the effective nucleon 
and dibaryon masses (3.7) and (3.8) and by the effective masses of the mesons (4.3) and 
(4.4). The graphical representations for the MF propagators are shown in Fig. 13. 

The total Green's functions can be determined self-consistently by solving a system of 
Gorkov-Dyson equations. As an example we derive an equation for the cr-meson propa- 
gator. Let us multiply Eq.(3.3) by o", take the time-ordered product of the equation, and 
find the average value of the equation over the ground state 

(-□ - ml) < Ta{l)a{2) >= 6\l, 2) + < T'l'(l)*(l)<7(2) > 
+2h^ < T{mD + K(j{l))^*{l)^{l)a{2) > . 
Taking into account the second-order terms with respect to the operator fields, we obtain 

(_□ _ ml) < ra(l)a(2) >= 5\1, 2) + 2hlp'i,s < Ta(l)a(2) > 
+ {2mhh^^s)i< mi)a{^) > + < T0*{l)a{2) >). 

The second term on the right hand side redefines the mass of the cr-meson (see Fig. 13 and 
Eq.(4.3)). In the momentum representation, Eq.(5.5) takes the form 



D--(A;) = D--(A;) + D''''{k)2m*j,K^pls{D''''{k) + D'^'^k)), (5.6) 

where 

D'^^ik) = — (5.7) 

is the (T-meson MF propagator. 

The equations for other Green's functions can be obtained in similar way to give 



D^^ik) = D'^'^{k)2m},h^^phs{D^,^{k) + Df-{k)), (5.8) 

D'^^ik) = b''''{k)2m*i,K^s{D'''''{k) + D^^{k)), (5.9) 
D'^'^-'ik) = D'^'^{k)2mlK^s{D''''''\k)+D^^\k)), (5.10) 



D^A^) = + b-^^{k)K^s[{'^^^D + k)rDr{k) 

+ i2,,h-k)rDf-{k)], 



(5.11) 
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Here 



+ {2fil-k)rD^'^'{k)], 

DW(fc) = D^^*{k) + Dw*(A;)[/i^^/;^(2/i^ + k^D^'P^k) 
+2mlh,^^D''^*{k)], 

D^'^'ik) = b^*^{k)[h^^sWD - k)rDr*ik) 
+2m*j,K^D'^^\k)]. 



D^^\k) 



_ ^2 
1 



5.12) 
5.13) 
5.14) 
5.15) 
5.16) 

5.17) 
5.18) 



are the cu-meson and dibaryon MF propagators. 

The system of equations (5.6) and (5.8)-(5.16) is shown graphically in Fig.l4. It admits 
an explicit solution. The propagators D'^'^* (k) and D'^*'^ {k) are expressible in terms of 
the propagators D'^'^* {k) and D'^'^* {k) and analogous propagators for the tu-meson. These 
propagators in turn are expressible through the propagators D"^"^* {k) and D'^*'^* (/c). 

The system of two equations for dibaryon Green's functions 



D'^^*{k) = b^^\k) + L)w*(A;)Sw*(A;)L)w*(A;) + L)w*(A;)sw(A;)L)^V*(A;), 
D^'^'{k) = b^''^{k)i:'^*^{k)D^*^*{k) + b^*^{k)i:^''^*{k)D^^*{k), 



(5.19) 



where 



S^^*(fc) = J:^*n-k) = {K^sYi'^^^D + k),D-::ik)i2fi}, + k), 

+ i2mhK^Yb'^-{k), 
S^^(A;) = S^>*(-A;) = {h^^sn'^f^h + k)^W^:{k){2ii*j, - k), 

+ {2mlh^^)'D-'^{k) 



(5.20) 



constitutes therefore a closed system of equations. These equations are shown graphically 
in Fig. 15. They are identical to the Gorkov equations in the theory of superconductiv- 



ity |^6[. The relativistic version of these equations for the SU(2) color quark matter is 
discussed in Ref. The system has solutions 

(Z)w*(A;)-i-Sw*(fc))(D^*<^(fc)-i-E^*^(A:))-S^*^*(A;)Sw(A;)' ^' ^ 
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/^From expression 



D^^(A;) = D^^*{k)i:^^*{k)D-^^*{k) + D^^*{k)E^^{k)D^*'^*{k) (5.23) 

we find also 

D^^ik) = . . (5.24) 

(Dw*(A;)-i - EW(A;))(L)vV(A;)-i - EW(fc)) - S^V*(A;)Sw(fc) 

The Green's functions for other particles are expressible in terms of the constructed 
dibaryon Green's functions. 

5.2 Dispersion laws for elementary excitations 

The mass operators J]'^'^* [k) has the form 



Denoting terms appearing in this expression, successively, by a, 6, c, and (i, one can 
write the following representations for the mass operators: S"^*^* [k) = a + b + c + d, 
S^>(A;) = a- b + c + d, and Sw(A;) = (A;) = a - c + ci. 

The chemical potential of the system is determined from the relation 

- m*i^ = S^^* (0) - S^'^(O) (5.26) 

which constitutes a relativistic extension of the non-relativistic relation derived by Hugen- 
holtz and Pines for Bose systems. For k = 0, b = c = 0, and the right hand side of 



Eq.(5.26) vanishes. Therefore, the MF relation (4.5) remains valid. A relation of such a 
kind is necessary to get a pole in the dibaryon Green's functions at a; = k = i.e. for the 
existence of sound in the medium. 

The spectrum of elementary excitations of the system is determined by zeros of the 
inverse Green's functions: 

{D'P'^* (k)-'^ - T.'P'P*{k)){b'P''P{kY^ - S'^*'^(A;)) - S'^*^*(A;)S'^^(A;) = 0. (5.27) 

To eliminate poles coming from the meson propagators, we multiply the denominator of 

- 2 

the Green's functions by {k"^ — m^){k'^ — m^)^. We then get an equivalent Ath order 
polynomial with respect to u 



2 . 



ELo^'"c„ = (P - ml)ie -rhD^x 
[(Z)w*(A;)-i _ Sw*(A;))(£)^V(A;)-i - S^*^(A;)) - S^>*(A;)Sw(A;)] 
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with coefficients 

Co = + k^iml + ml) 

+8k\rhli,*^hlp'hs - rhlfi*^hlp'hs), 
ci = -4k6 - k4(3m2 + 37hl + ) 
-2k^{mlml + 2ml,j,*^ + 2mli^*^ + An*j^hlp'i,s - ^f^Dhlp^s) 39) 
-ArfilmliJ,*^ + 8rfilii*^hlp'hs, 
C2 = 6k4 + k^{3rhl + 3ml + ^I^d) 
+mlml + Amln*^ + Amlii*^ - 8ii*^hlp%s, 
C3 = -4k2 -ml -ml- A/i*^, 
— 1. 

The polynomial (5.28) determines four different excitations, two of dibaryon type (par- 
ticles and antiparticles, the first excitation is sound) and two of a- and a;-meson types. 
These four poles occur in all Green's functions because of the a-uj-(p-(p* mixing. Such a 
mixing occurs because the a- and ci;-mesons can be absorbed by the dibaryons in the con- 
densate as a result of which the dibaryons leave the condensate and propagate as normal 
particles. The mixing describes also processes with creation of dibaryon- antidibaryon pairs 
with subsequent absorption of dibaryons by condensate and propagation of antidibaryons. 

One can verify that the Green's function D"'"'{k) has no pole at k'^ — ml. The poles 
of D'^'^{k) coincide with the poles of the dibaryon Green's functions. 

The second term in the a;-meson propagator in Eqs.(5.11) is factorizable and there 
is no Qfjtv tensor structure. Therefore, the pole at k"^ — in the first term cannot be 
cancelled by the second term, as in case of the cr-meson. The a;-meson in heterophase 
nucleon-dibaryon matter has therefore two branches of excitations. 

The velocity of sound can be found from Eq.(5.27) by keeping terms of order 0{\s?) 
and 0(0;^). In this hmit only the sound mode to — Ugik) survives. We get 

= ( ) k=o = (5-30) 



dk 

with 



l + a 



a^2p%s'^{^-^). (5.31) 

When the condition for stability (3.21) is fulfilled, the value is positive and less than 
unity (i.e. less than the velocity of light). 

The numerical analysis of Eq.(5.27) shows that there are no complex or negative uj'^ 
when the inequality (3.21) is satisfied. This result implies the stability of the ground state 
of the system for al > 0. If the poles occur symmetrically at ± u!a{k), the denominator of 
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the dibaryon Green's function has no zeros on the imaginary axis of the complex cj-plane. 
One can check that after a Wick rotation, the denominator really becomes positive for all 
values of uj and k. 

The group velocities of all four excitations are less than the velocity of light in absolute 
values. This is quite natural, because the model is relativist ically invariant. It the limit 
k — >■ oo we can keep the leading terms in k^ in Eq.(5.28). The polynomial (5.29) can then 
be summed up to (cu^ — k^)^. The dispersion laws for all four types of excitations behave 
asymptotically as u;^ ~ k^, and the group velocities approach unity in the limit of large 

k2. 



6 DISCUSSIONS AND CONCLUSIONS 

The qualitative estimates based on a model of non-interacting nucleons and dibaryons 
show that in normal nuclear matter a dibaryon Bose condensate does not exist provided 
the inequality m^) > 1.96 GeV is fulfilled. A more accurate estimate can be made on 
the basis of the relativistic MF model. /^From the requirement of absence of a dibaryon 
Bose condensate for pxv < Po? where po = 0.15 fm~^ is the saturation density for nuclear 
matter, we get for /i^ = 2(7^^ a constraint 

mr, > 1.89 GeV. 

This constraint is valid provided that dibaryon matter is stable against compression. It 
follows that the rfi-resonance with a mass rriD = 1.92 GeV does not affect the properties 
of ordinary nuclei. 

It would be interesting to check astrophysical data for the presence of a dibaryon 
condensate in the interiors of massive neutron stars as well as possible signatures of 
their instability caused by dibaryons. From the existence of massive neutron stars, one 
can put a lower limit on the masses of dibaryons. The estimates ( p. 21) and the results 
based on a more realistic model including interactions between dibaryons show that 



such a constraint can be physically significant. Conversely, the experimental discovery of 
dibaryons will have important astrophysical implications. 

Phase transitions of nuclear matter to strange matter |^5|, |5^ have been widely dis- 



cussed in the literature (for a review see fS^). Dense nuclear matter with a dibaryon 



Bose condensate can exist as an intermediate state below the quark-gluon phase transi- 
tion. This is the case when dibaryon matter is stable against compression. If dibaryon 
matter is unstable against compression, the creation of dibaryons could be a possible 
mechanism for the phase transition to quark matter. 
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The soft core of the ifif-interaction ||3^ is responsible for the relatively low value of 
the critical density for formation of iJ-dibaryons in nuclear and neutron matter and for the 
possible instability of if-matter. The energetically favorable compression of if-matter will 
eventually lead to the formation of absolutely stable strange matter. Possible astrophysical 
examples are bursters and roentgen pulsars which accret matter from companion stars. 
This leads to an increase of the mass and central density of these neutron stars. Once the 
density exceeds the critical density, if-particles can be created, leading to the formation 
of strange matter. The neutron star converts then to a strange star ||5^, |59|, Q . 

The experimental observation of a Bose condensate of dibaryons in heavy-ion collisions 
would be of great importance for understanding physics of nuclear matter at supranuclear 
densities. In the center-of-mass frame of the condensate a large fraction of dibaryons 
has zero velocities. When the density decreases, dibaryons in the condensate decay to 
their specific channels. Experimentalists would observe in every collision events with the 
same invariant mass rriD and the same total momentum. An excess of such events can be 
considered as a possible signature for the formation of a dibaryon condensate in heavy-ion 
collisions. 

The contribution of dibaryons out of the condensate to pressure and energy density 
can be calculated in relativistic Hartree approximation only. We constructed Green's 
functions of the system. The one-loop calculation of the EOS for heterophase nucleon- 
dibaryon matter will be given elsewhere 
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FIGURE CAPTIONS 



Fig.l. A schematic representation for the process of dibaryon Bose condensation in 
neutron matter in the ideal gas approximation. When the chemical potential of nucle- 
ons exceeds mo/'^-i production of dibaryons becomes energetically favorable. Dibaryons 
are Bose particles. They are accumulated in the ground state of the system with zero 
momentum and form a Bose condensate. 

Fig.2. Schematic behavior of the pressure as a function of the density in the ideal gas 
approximation. Above the critical density for formation of the dibaryon Bose condensate 
(indicated by the arrow), the neutron chemical potential /x„ is frozen at a value m£)/2 and 
the neutron density is fixed, but the total baryon number density is still increasing along 
the dashed line due to the dibaryon formation. Dibaryons are in the ground state and 
do not contribute to the pressure. Therefore the pressure for the binary mixture remains 
constant (horizontal dotted line). The solid line gives the EOS for homophase neutron 
matter. 

Fig. 3. Diagrams contributing to the coupling constants of the oo - and o"-mesons with 
dibaryons in the additive model. The dibaryon D couples strongly to the iViV-channel (a), 
the (i'-dibaryon decays into the TriVA^-channel {h). The a;-meson is coupled to nucleons 
(a, 6), the cr-meson is coupled to nucleons and the 7r-meson in the d' {h). 

Fig.4. The a- and 7r-meson exchange current contributions to the cr-meson coupling 
constants with the dibaryon D coupled strongly to the A^A^-channel (a), and with the 
d'-dibaryon decaying into the TrA^A^-channel (6). 

Fig.5 (a,b). The critical density for occurrence of a Bose condensate of dibaryons in 
nuclear (a) and neutron (6) matter versus the cr-meson coupling constant /ig. for mo = 1.88 
GeV (= 2miv; the long dashed curve No. 1), 1.96 GeV (the sohd curve No. 2), etc. with a 
step 80 MeV. The results for the di(1920) and d'(2060) dibaryons are shown (the dashed 
curves). Dibaryon matter is stable against compression when the square of the sound 
velocity is positive. This is the case for hcr/{2g^) < 0.8754. The value po — 0.15 fm~^ 
is the saturation density for nuclear matter. The occurrence of if-dibaryons in nuclear 
and neutron matter is denoted by the crosses. 

Fig.6. Possibilities for the behavior of the difference 2//jv— A*d between the two-nucleon 
and the dibaryon chemical potential versus the dibaryon fraction 2p'jjy/pTv- Figure (a) 
shows unstable equilibrium states. Figure (6) shows stable equihbrium states. 

Fig. 7 (a,b) The effective nucleon mass in GeV versus the dibaryon fraction 
'^Pdv/Ptv in heterophasc matter (a). The results do not depend on the dibaryon mass. 
The difference 2/i7v — Pd between the two nucleon chemical potentials and the dibaryon 
chemical potential versus the dibaryon fraction 2p'j^y/pTv (&)• The results are given for 
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total baryon densities 1, 2, 3, 4, 5, and 6 times the saturation density po- The normal 
homophase matter is stable when 2/1^ — < and p'fjy = 0. An intersection of a curve 
with a negative slope with the horizontal line 2^j^ — fi^ = indicates occurrence of a 
stable equilibrium in hcterophase matter. Two such states for nuclear and neutron matter, 
occurring at Ptv = 3po and Ptv = 2poj are denoted by the arrows. The results are given 
for nir) = 1.96 GeV. The dibaryon mass does not enter the self-consistency condition 
(4.11) and enters linearly in the difference 2/ijv — Pd, and the curves for other dibaryon 
masses can be obtained by vertical parallel displacements. The results for m£, = 2.06 GeV 
((i'-dibaryon) can be obtained e.g. by a 100 McV negative shift, etc. The solid lines stand 
for nuclear (7 = 4) matter, the dashed lines stand for neutron (7 = 2) matter. 

Fig. 8. The energy per nucleon in homophase nuclear and neutron matter (solid lines) 
and in heterophase matter (dashed and dotted lines) versus the total baryon number 
density pTv = Pnv + '^Pdv^^^ (ii(1920) and ci'(2060) dibaryons using hu, — 2g'^and 
ha/{2ga) — 0.8. The dibaryon Bose condensation decreases the energy of the ground 
states. The occurrence of i?-particles in nuclear and neutron matter is shown. It results 
in the formation of strange matter. The value po is the saturation density of nuclear 
matter. 

Fig. 9. Equation of state for homophase nuclear and neutron matter (sohd lines) and 
for heterophase nuclear- and neutron-dibaryon matter (dashed and dotted lines) versus 
the total baryon number density for di(1920) and d'(2060) dibaryons at — 2g'^and 
ha/ i'^Qa) — 0.8. Dibaryon Bose condensation at high densities softens the EOS for nuclear 
and neutron matter. The critical density for the occurrence of if-particles in nuclear and 
neutron matter is indicated. 

Fig.lO. There are two kinds of vertices corresponding to interactions of the lo- and 
(7-mesons with dibaryons and a vertex describing creation and absorption of dibaryons by 
the condensate. 

Fig.ll. Pictorial representation of the series for the uj- and ci-mcsons mean fields 
(Eqs.(4.1) and (4.2) ). The diagrams can be summed up to produce (i)the dressed meson 
MP propagators without modification of the meson vertices or {ii)the dressed meson 
MP vertices with dibaryons without modification of the meson MP propagators and the 
meson-nucleon MP vertices. 

Fig. 12. The dressed uj- and cr- meson MP vertices with dibaryons. 

Fig. 13. The Dyson equations for the MP propagators of the a- and cj-mcsons, nu- 
cleons, and dibaryons. Thin lines define the bare propagators, thick lines define the MP 
propagators. The dashed lines, the solid lines, and the double solid lines describe, respec- 
tively, the meson propagators, the nucleon propagator, and the dibaryon propagator. 
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Fig. 14. Pictorial representation of the Gorkov- Dyson equations for the complete 
Green's functions in the heterophase nucleon-dibaryon matter. Note the correspon- 
dence between the diagrams and the equations in the text: (a)- Eqs.(5.6) and (5.11), 
(6)- Eq.(5.8), (c)- Eqs.(5.9) and (5.12), {d)- Eqs.(5.10) and (5.13), (e)- Eq.(5.14), (/)- 
Eq.(5.16), and (^)-Eq.(5.15). 

Fig. 15. Pictorial representation of the system of two coupled equations (5.19) for the 
normal and anomalous dibaryon Green's functions. 
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